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ABSTRACT 

The nearest accretion disc to us in space if not time was the protosolar nebula. 
Remnants of this nebula thus potentially offer unique insight into how discs work. In 
particular the existence of chondrules, which must have formed in the disc as small 
molten droplets, requires strong and intermittent heating of disc material. We argue 
that this places important constraints on the way gravitational energy is released in 
accretion discs, which are not met by current shearing-box simulations of MRI-driven 
dynamos. A deeper understanding of accretion energy release in discs may require a 
better model for these dynamos. 
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1 INTRODUCTION 

It is now generally agreed that accretion discs are driven 
, mainly by magnetic torques (Shakura & Sunyaev, 1973) and 
■ that the magnetic fields in such discs are maintained by 

local dynamo processes. What is not understood however, 
, despite considerable theoretical efforts, is how such dynamo 
' processes work, and exactly where and in what form the 
' accretion energy is released (King et al., 2007; Blackman, 

2010). 

Given the difficulty of this problem, we should look for 
. observational evidence which bears on it. We have unique 
' insight into the remains of one accretion disc - the protosolar 
nebula. It is therefore worth asking if the present-day solar 
system offers evidence which can constrain the way energy 
is released in accretion discs. 

Among the oldest solid objects in the solar system are 
chondrules, the round grains present in the majority of mete- 
orites. These must have formed as molten droplets in space 
before being accreted. The need to heat them sufficiently 
implies a connection with energy release in the protosolar 
nebula. We consider here the overall energetics required of 
the heating process. We argue that since most of the mete- 
oritic material has been subject to this heating, the energy 
source required must be quite widespread and substantial, 
and must therefore involve the major local source of energy, 
i.e. disc accretion. We find that around 10 per cent of the 
accretion energy is required. If so, it is evident that the exis- 
tence of chondrules provides fundamental information about 
energy release in accretion discs. 

This paper is organised as follows. In Section 2 we give 
a brief introduction to chondrules and ideas about their for- 



mation. In Section 3 we consider the global energetics of 
the heating process. A likely mechanism for the provision 
of transient heating within the solar nebula involves strong 
shocks (Deschl & Connolly, 2002, Connolly et al., 2006). In 
Section 4, we briefly present the model of chondrule forma- 
tion through shock heating suggested by Deschl & Connolly 
(2002) but argue that their proposed mechanism for gener- 
ating these shocks by gravitational instabilities within the 
disc requires disc properties which do not sit easily with our 
current understanding of disc evolution. In Section 5 we de- 
scribe the disc properties which appear likely to hold at the 
time of chondrule formation. The need for significant disc 
dissipation in the form of shocks, together with evidence for 
magnetic fields within the nebula at that time (Section 6) 
lead us to a picture of a disc dynamo (Section 7). This differs 
markedly from the results of current numerical simulations. 
Our picture draws on the dynamo model of Baggaley et 
al. (2009a, b) and involves thin flux ropes and reconnection, 
analogous to models of flaring on the solar surface. Section 8 
is a discussion. 



2 CHONDRULES 

Most meteorites are chondrites (more than 75 per cent. Sears 
& Dodd, 1988; Hutchinson 2004), and the most abundant 
constituent of the majority of chondritic meteorite groups 
are chondrules (Grossman et al., 1988). Chondrules are small 
particles of silicate material (typically around a millimetre 
in size, Weisberg et al., 2006) that experienced melting be- 
fore incorporation into chondritic meteorite parent bodies 
(Hewins, 1996). The properties of chondrules are consistent 
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with formation in the protosolar nebula from primitive ma- 
terial, and thus provide an important record of events in this 
disc. Isotopic dating of chondrules indicates ages of 4.5 Gyr, 
and also suggests that high temperature nebular processes, 
such as chondrulc formation, lasted for about 3 — 5 Myr 
(Russell et al., 2006; Scott, 2007). This timescale coincides 
with the observed lifetime of protostellar discs around solar- 
type stars (e.g. Hernandez et al., 2009). Chondrules formed 
a few Myr after the earliest CAIs (Ca-Al-rich inclusions) 
(Swindle ct al., 1996). Those formed at lower temperatures 
(~ 1500 K) and in general have longer cooling times (Davis 
& MacPherson, 1996), and are part of the same sequence of 
heating events. 

The source of most meteorites is now recognised as the 
asteroid belt (Hutchinson, 2004) . The present asteroid belt is 
thought to consist of a tiny fraction of the planetesimals that 
were originally at these radii, the majority being expelled 
by dynamical resonances once the large planets (Jupiter 
and Saturn) had acquired sufficient mass (Morbidelli et al., 
2005). Thus we may regard the chondritic meteorites as a 
trace sample of what remained in the protoplanetary disc at 
radii around 3 AU once disc evolution had finished, much 
of the solid material had been accreted into planets, and 
any gas not accreted to form the giant planets had been 
expelled. Unless the dynamical processes which expelled the 
bulk of the solid material from the region of the asteroid belt 
had some means of preferentially leaving behind chondritic 
material, it follows that the chondritic meteorites give some 
idea of the planetesimals formed in that region, and there- 
fore of the physical processes in the disc towards the end of 
the major accretion phase. Thus, although there are models 
of chondrule formation which involve recycling of material 
radially in the disc, either by diffusive processes (e.g. Mor- 
fiU 1983, Clarke & Pringle, 1988) or, more plausibly, by a 
more direct route (Shu et al., 2001; see also King, Pringle & 
Livio, 2007), we shall assume here that chondrule formation 
took place more or less tn sttu and ask what implications 
this might have for the workings of an accretion disc. 

Many different mechanisms have been put forward for 
the formation of chondrules (see, for example, the review 
by Boss, 1996). As discussed there, the main problem with 
the formation of chondrules in situ is that the temperatures 
needed in their formation (1500 - 2000 K, Hewins, 1988; 
Hutchinson, 2004; Connolly et al., 2006) greatly exceed the 
likely ambient temperatures in the protoplanetary disc at 
late times. Higher disc temperatures are possible at earlier 
times when the disc was more massive, disc evolution more 
rapid, and the accretion rates large, but chondrules formed 
at that time would mostly have been carried inwards by 
the accretion disc and deposited into the Sun. Chondrule 
formation in a protostellar disc as it reaches its later stages of 
evolution requires a series of transient heating events within 
the disc (Connolly et al., 2006). Thus 

'Although a very large number of mechanisms have been 
proposed for making chondrules, only one has been shown 
to satisfy most of the constraints: shock heating of nebula 
gas and entrained silicates.' (Scott, 2007). 

However the same author also remarks that 

'There is no observational evidence for shocks in pro- 
toplanetary disks, and a single source of shocks that could 
have operated over several million years has not been demon- 
strated.' (Scott, 2007). 



We attempt this in the present paper. 



3 ENERGETICS OF THE FORMATION 
PROCESS 

As remarked above, chondritic meteorites appear to come 

from the asteroid belt. Accordingly we assume that in situ 
chondrulc formation occurs in a region of the protostellar 
disc at radii in the range 2-5 AU, and so at around a 
typical radius of 

i? = 3AU = 4.5 X lO^^cm. (1) 
At this radius disc material has circular velocity 

and thus a specific kinetic energy 

i.„..iv^.I.5xI0"(-M.)(3|5)e,gg-. (3) 

We can think of this as being equivalent to the accretion 
energy available locally. This is the energy source that we 
shall need to use to power chondrule formation. 

The specific energy required to heat the material from 
a few hundred K up to ~ 1800 K and then melt it is given 
by Wasson (1996) as 

S^eit = 2.1 X 10'°ergg-\ (4) 

To form chondrules it is necessary to heat the proto- 
chondrules (or 'dust balls') to melting point. The fraction of 
chondrule to matrix (the remainder of the chondrite, which 
has never melted) in chondrites is laxge, and in addition 
there is evidence that about 15 per cent of chondrules have 
undergone more than one molting event. Thus wo shall as- 
sume that on average each proto- chondrule (or 'dust-ball') 
is subject to around one melting event. 

Since _Ekin ^ E^^it, it seems there is plenty of energy 
available. However, this is not the whole story. Once heated 
to melting point, the chondrules must cool. If they simply 
radiated their energy into space as black bodies they would 
cool on a timescale 

~ Wrf^' 

where a chondrule is assumed to have mass m, radius r and 
temperature T, and a is the Stefan-Boltzmann constant. 
Thus a typical cooling time for a chondrule would be (see, 
for example, the cooling curves given by Wasson, 1996) 

^'^-^ = (2.1xl0i0erg/g) (l^) (ij/^) (i800k) 

where pchon is the density of the chondrule. 

In practice, in order for the chondrules to have the prop- 
erties observed, they must cool on a much longer timescale 
than this, i.e. tcooi >■ trad (Scott 2007). Estimates of cool- 
ing times vary greatly. Levi & Araki (1989) suggest tcooi ~ 
2 X 10^ to 2 X 10" s. Hutchinson (2004) agrees with this at 
the short end, but indicates that some cooling times must be 
much longer, in the range tcooi ~ 2 x 10^ to 3 x 10® s. From 
the detailed discussion given by Desch & Connolly (2002) it 
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is evident that the cooling histories of chondrules are var- 
ied and complicated. Scott et al. (1996) argue that the disc 
must have been 'a maelstrom where temperatures fluctuated 
through ~ 1000 K and solids experienced "multiple cycles 
of molting, evaporation, recondensation, crystallization and 
aggregation" '. 

Desch & Connolly (2002) give characteristic cooling 
times in the range of hours to days. This is in line with 
the recent review by Connolly et al. (2006), who quote the 
majority of chondulos as having a cooling time of ~ 10^ s. 
Thus wo adopt a typical cooling time for chondrules of 



icool — 10' s. 



(7) 



The fact that tcooi ^ trad has implications for the en- 
ergetics of the formation process. Not only must the chon- 
drules be heated to melting temperature, but they must also 
be kept near that temperature for a prolonged period. This 
implies that not just the proto-chondrulcs but also the ma- 
terial surrounding them must be heated to around 2000 K 
(Wasson, 1996). In addition, a sufficient volume of surround- 
ing material must be heated so that its optical thickness im- 
plies a cooling time of hours to days. The specific energy 
required to heat the surrounding gas to a temperature of 
T = 2000 K is given approximately by 



2mH 



^ (soSok) 



(8) 



where we have assumed that the typical gas particle is a 
hydrogen molecule, and have ignored molecular dissociation. 
This implies, again assuming that on average each element 
of chondritic material is subject to one melting event, that 
a fraction / of around 

/ ^ « 0.08 (9) 

-C'kin 

of the available accretion energy must be used to provide 
the heating process for chondrule formation. 



4 SHOCK HEATING 

Desch & Connolly (2002) have proposed a detailed model for 
the formation of chondrules in terms of the thermal process- 
ing of particles in shocks in the accretion disc. Their model 
takes account of the properties of the shock, including gas- 
drag heating of particles which are not slowed instantly in 
the gas shock, radiative processes including the disassoci- 
ation and recombination of H2, and particle evaporation. 
Their model can account for the thermal history of chon- 
drules. Their canonical model requires a shock velocity of 



7 X 10^cms"\ 



(10) 



For their assumed ambient disc temperature of T = 300 K 
and so sound speed (assuming a gas of molecular hydrogen) 
of Cs = 1.5 X 10^ cm s^^, this corresponds to a Mach number 
of around ~ 5. 

However, while this model can successfully account for 
chondrule formation in the disc, the proposed source of the 
shocks within the disc does not sit easily with our current 
understanding of protostellar disc evolution. Desch & Con- 
nolly (2002) propose that the shocks are caused by grav- 
itational instabilities within the disc. But this proposal is 



at odds with their disc properties, which are taken from a 
disc model by Bell et al. (1997) with an accretion rate of 
M = 10~*Mo yr~^ and a viscous parameter of a = 10""'. 
The disc density is assumed to be p = 10~® g cm~^, which 
for a disc semi-thickness H = (cs/V^)-R ~ 4 x lO"'^ cm gives 
a disc surface density of S = 8000 g cm~^. This gives rise to 
a strongly unstable disc with a Toomre parameter (Toomre 
1964) of 



Q 



ttGS 



: 3.4 > 1. 



(11) 



Such strong gravitational instability is indeed required if this 
process is to give rise to shocks with velocities Vs ~ O.4V0. It 
is, however, hard to envisage how all the material in such an 
strongly unstable disc would manage to avoid being strongly 
shocked a large number of times. 

Such a disc has a mass at radii around 3 AU (say, at 
radii 2-4 AU) of 



Mdisc(3AU) w nR H w O.O25M0, 



(12) 



Such a violently unstable disc is likely either to fragment 

into predominantly gaseous bodies, gathering any chon- 
drules with them, or, if not, give rise to strong gravita- 
tional torquing and so a high accretion rate. The accretion 
rate expected for a disc which is so gravitationally unsta- 
ble is the same as one for which the viscous parameter has 
a « 0.06 (Lodato & Rico, 2004; Rice, Lodato & Armitage, 
2005) which would give (Pringle 1981) 



M = « 4 X lO"^M0yr"\ 



(13) 



Such a high accretion rate would give a midplane disc 

temperature of ~ 2000 K (Bell et al., 1997) and contradicts 
the original assumption that M = W~^Mq yr^^. It would 
give a local disc lifetime of around Mdisc/M « 600 yr. 

Similar considerations of accretion rates and lifetimes 
apply to the more detailed models presented by Boss & 
Durison (2005; see also Boss, 2007). These authors stress 
that early times, when infall onto the disc was high, the disc 
is most likely massive enough for self-gravity to play a ma- 
jor role in angular momentum redistribution (see, for exam- 
ple, Lin & Pringle, 1990). They also note, however, that the 
local, non-axisymmetric instabilities driven by self-gravity 
occur at, or around, co-rotation (see also the discussion in 
Cossins et al., 2009). This implies that in general, through- 
out the bulk of the disc, the relative velocity between disc 
gas and the (unstable, but transient) spiral pattern is small. 
Thus, in general, the short-lived, localised, transient spirals 
which are the typical outcomes of self-gravitational disc in- 
stabilities, are not able to drive strong shocks in the gas. In 
this context, in the numerical simulations of Boss & Durisen 
(2005) it is only close to (within a few grid cells of) the inner 
grid boundary that they are able to achieve a shock pattern 
which is neither close to co-rotation, nor strongly trailing. 

Thus while we cannot rule out the possibility that a 
massive disc and self-gravitational instabilities arc able to 
provide the necessary shocks, we now consider what we re- 
gard as more likely properties for the disc in the region of a 
few AU at the time of chondrule formation. 
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5 DISC PROPERTIES 

From timings deduced from isotopic data, Russell et al. 
(2006) conclude that chondrule formation lasted for around 
3-5 Myr. This is also the timescale on which young 
solar-type stars appear to lose their protoplanetary discs 
(Hernandez et al. 2009). Thus it seems reasonable to con- 
clude that chondrules were formed during the final stages of 
the protoplanetary disc when the accretion rates were low 
(M < 10"'^ Mq yr-\ Ciesla & Charnley, 2006). 

Models of protostellar discs at low steady accretion 
rates arc given by Bell ct al (1997). At the accretion rates 
we consider, it is likely that heating from the central star 
provides a non-negligible, and perhaps dominant, contribu- 
tion in the surface layers, which may cause some changes 
on the midplane (e.g. Chiang et al., 2001). There is also 
the possibility that such discs contain 'dead zones' (Balbus 

6 Hawley 2000; Terquem, 2008) where the MRI is unable 
to operate because of high magnetic diffusivity caused by 
low ionization. However for the canonical disc parameters 
we use (see below) Matsumura & Pudritz (2003) argue that 
the MRI should be fully operational. This is in line with 
our argument above that the full accretion energy must be 
accessed to power chondrule formation. 

At the time of chondrule formation wc take the disc 
mass as a few times the minimum mass solar nebula (say, 
~ O.O6M0) and its age as ~ 3 Myr. An accretion rate of 
order M ~ 3 x 10~®Mq yr~^ is appropriate. Given this, the 
models of Bell et al (1997) suggest the following parameters 
for the disc in the region of ii « 3 AU. The disc (midplane) 
density is 



Pdisc 



10 gem , 



and the midplane temperature is 

Tdisc = 200K, 

with disc thickness 

H = 0.05R = 2.2 X lO^^cm 

and viscosity parameter 

a = 0.01. 



(14) 



(15 



(16) 



(17) 



In this model most of the disc mass still resides at radii 
10 - 20 AU, and it is these radii which control the evolution 
timescale of the disc, and so the accretion rate. These pa- 
rameters arc in line with the fully ionized models of Terquem 
(2008) and with the disc properties assumed by Ida & Lin 
(2004) in their modelling of planet core formation. For these 
parameters the disc sound speed is 

c, = 1.2 X 10'''cms"\ (18) 
and the gas pressure is 

= 0.83 dyne cm"^ (19) 



6 MAGNETIC FIELDS 

Remnant magnetic fields found in chondritic material imply 
cooling through the Curie temperature of ^ 600 K in the 
presence of a magnetic field. This offers clear evidence that 
chondrule formation took place in regions with magnetic 
fields in the range B ~ 1 - 10 G (Levy & Sonett, 1978; Levy 



& Araki, 1989) . This is close to the maximum (cquipartition) 
field strength B of a magnetic flux tube in an accretion disc 
- equating the magnetic pressure /8n with the local gas 
pressure pg for the disc parameters given above gives 



4.6 G. 



(20) 



We note that to achieve the strength of viscosity {a = 0.01) 

assumed above, wc only require an average magnetic field 
(B) within the disc given by (Shakura & Sunyaev, 1973) 



8-!TPg 



1/2 



(21) 



(22) 



This suggests that chondrule formation could have taken 

place in a disc in which a dynamo mechanism driven by the 
magneto-rotational instability was operating. 



7 THE DYNAMO PROCESS AND 
CHONDRULE FORMATION 

We have seen that chondrule formation requires the heating 
of each element of the disc (gas plus proto-chondrules) up 
to a temperature of around T « 2000 K « lOTdisc, for a 
time ~ tcooi ~ 10'"' s, on average once. This must happen 
during the time that this material is at a radius of order 3 
AU. Once the gas is at this temperature it becomes over- 
pressured relative to its surroundings, and so will tend to 
expand, and therefore cool adiabatically, at its now elevated 
sound speed of c* ~ 3.8 x 10^ cm s~^. For this adiabatic 
cooling to take longer than tcooi we require the size of the 
heated region h to be such that 



h ^ Cg X ^cool 



! 3.8 X 10^°cm. 



(23) 



We have also noted (Deschl & Connolly, 2002) that a 
way of achieving such heating is through shocks within the 
disc, with shock velocities Vs ^ 7 x 10^ cm s"'^ ~ 6cs. 
However, although the standard shearing-box models of 
MRI-driven dynamos in accretion discs (see Davis, Stone & 
Pessah, 2009) can produce values of a of order 0.01, they 
predict velocities within the disc which are for the most 
part strongly subsonic. In this picture, chondrule formation 
would thus remain unexplained. 

However, there are indications that such models do not 
give a complete description of dynamo processes in accretion 
discs (cf King, Pringle & Livio, 2007). One problem is that 
these models use the MHD approximation with a standard 
(isotropic) magnetic difi^usivity 77 and standard (isotropic) 
Navier-Stokes viscosity v. Fromang et al. (2007) demon- 
strate that the results obtained depend on the assumed 
numerical parameters: the Reynolds number TZ = CsH/v, 
the magnetic Reynolds number TZm = CsHrj and, criti- 
cally, on the ratio of the two, the magnetic Prandtl num- 
ber Vm = TZ-m/TZ = vlr\. There are also other physical 
processes which can provide complications such as the Hall 
effect (Balbus & Hawley, 2000), and anisotropic transport 
processes (e.g. Dong & Stone, 2009). Schekochihin and co- 
workers (e.g. Schekochihin et al., 2004, 2005) have also ar- 
gued that the structure of a magnetic field in a turbulent 
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flow at such high Reynolds numbers can depend criticaUy 
on the Prandtl number, and that this has imphcations for 
dynamos under such conditions. In addition, Heitsch et al. 
(2008) and Zweibel & Heitsch (2008) have discussed the im- 
phcations for magnetic growth and structure in turbulent 
media where ambipolar diffusion plays a significant role. 

7.1 An alternative dynamo picture 

In the light of the above, we speculate how an accretion disc 
dynamo might produce the intermittent high-energy events 
which appear to be demanded by chondrule formation. As 
an example we consider a new model for a dynamo in a 
turbulent medium suggested by Baggaley et al. (2009a,b). 
The extent to which such a model is applicable to cooler 
accretion discs, such as the solar nebula, is somewhat un- 
certain, and we introduce it here because it is able to illus- 
trate the kind of properties with regard to energy release 
that we are looking for to facilitate chondrule formation. 
In this model the magnetic field is mainly confined to thin 
flux ropes which are advected by the flow. Magnetic dissi- 
pation only occurs via reconnections of the flux ropes, and 
so the magnetic dissipation is highly localised. This model 
can be viewed as an implementation of the limiting regime 
of infinitely large magnetic Reynolds number: magnetic dis- 
sipation can be safely neglected at all scales, but plays a 
crucial role through reconnection of field lines in permitting 
rearrangement of the field topology. Such rearrangements of 
the field topology result in conversion of magnetic energy 
into kinetic energy of the fluid, and thence to dissipation as 
heat. This contrasts to the usual models in which magnetic 
diffusivity converts magnetic energy to heat directly. 

With this in mind, the picture we propose is that the 
magnetic fleld in the disc can be thought of as a collection 
of loops of flux ropes. These loops are continually stretched 
by the azimuthal shear flow. The stretching increases the 
magnetic energy associated with the loop and at constant 
total (gas plus magnetic) pressure decreases the mass den- 
sity along the fleld line, and so increases the Alfven speed 
along the loop. To maintain an equilibrium distribution of 
loop properties, the loops must also continually undergo re- 
connection events. The crucial property of these events is 
that while each reconnection event in itself releases a neg- 
ligible amount of energy, it does release the magnetic fleld, 
which is then able to reconflgure itself (at the Alfven speed) 
and in so doing, heat the gas. The picture here thus differs 
fundamentally from that presented by Sonett (1979) and 
Levy & Araki (1989). They considered magnetic reconnec- 
tion in low density regions far from the disc plane (in the 
disc corona) and took only the magnetic energy dissipated 
by the reconnection events into account. We are assuming 
here that most of the accretion energy is released in the 
bulk of the disc. We should note, however, that some au- 
thors have suggested that a substantial fraction of the ac- 
cretion energy might be released in such low density regions 
(Tout & Pringle, 1990; Uzdensky & Goodman, 2008), ar- 
guing that magnetic buoyancy can advect energy efficiently 
away from the disc plane. For our picture to work, however, 
we require that a significant amount of reconnection occurs 
close to the disc plane, so that the bulk of the energy release 
occurs there and disc gas, along with the chondrule, can be 
efficiently shock heated. 



What we therefore require for chondrule formation is 
that some reconnection events occur in regions with suffi- 
ciently high Alfven speeds that shock velocities of Va « 6cs 
can be generated. For this to occur we require that a mi- 
nority of flux tubes have field strengths of order Boq and 
loop mass densities around 40 times lower than the mean 
disc density before reconnection occurs. Unfortunately the 
'fluctuation dynamo' model of Baggaley et al (2009a, b) is 
currently computed only in an incompressible medium. Fur- 
ther consideration of the model will be required before it is 
possible to establish whether or not such high-energy recon- 
nection events are likely in an accretion disc. 



7.2 Implications for magnetic field structure and 
internal disc properties 

In this picture we suppose that from time to time a recon- 
nection event occurs which causes a sufficiently large and 
rapid adjustment of fleld topology that a region of disc gas 
of size h is subject to shock heating, with velocity Vs « 7 km 
s~^. We need this heated region to stay sufficiently hot for 
a time tcooi ~ 10^ s. We can use this information to deduce 
requirements for the properties of the magnetic loops. 



7.2.1 Size of heated regions 

We expect the cooling timescale r for a region of size h and 
temperature T to be given by 



heat content 

T ~ . 

heat loss rate 
We expect that the heat content is simply 

heat content oc pTh^ , 



(24) 
(25) 



where p is the gas density. If the region is optically thick, so 
that heat loss is mainly by radiative transfer, then we expect 



heat loss rate cx x h , 

pKh 

where k is the opacity. Thus we expect 

2 7,2 

p Kh 



rp3 



(26) 



(27) 



For an accretion disc in thermal equilibrium, the rel- 
evant lengthscale is the disc thickness H, and the cooling 
timescale is the thermal timescale tth, given by (Pringle 
1981) 



tth ~ ~ 2.6 X 10' 
ail 



Vo.oiy 



(28) 



Using this we can obtain a very rough estimate of the 
size required for the shock-heated regions. To keep things 
simple we assume that the opacity k ~ const, (e.g. Bell et 
al., 1997), although if some chondritic material is vaporised 
and /or H2 is significantly dissociated, this might not be the 
case. Then from the computations of Deschl & Connolly 
(2002) we note that once the gas has relaxed behind the 
shock, the temperature is T ~ 1000 K, and so a factor of 
~ 50 above its pre-shock value, and similarly that the den- 
sity p has increased by a factor of order ~ 10. Thus using 
Equation [27] we find 



6 A. R. King and J. E. Pringle 



^^'-^x {50r X (10-^), (29) 
which gives 

h ~ 0.2H ~ 4 X lO" cm. (30) 



7.2.2 Heating events and loops 

We have seen that a single heating event is likely to have 
a size of around /i ~ 4 x 10^^ cm. Evidently this must be 
roughly the size of the reconnecting loops. This means that 
the rate per unit volume at which energy is dissipated locally 
by the shock is approximately 



gshock ~ X 



(^) 



(31) 



We can contrast this with an accretion disc in local 
thermal equilibrium, where the much lower average volume 
rate at which energy is dissipated within the disc is approx- 
imately 



tth 



(32) 



We can now use this to estimate the probability p that 
a random disc element is hot at any one time. The average 
rate at which shocks heat the disc is 



{qs) ^ p X gshock. 



(33) 



We have argued earlier (eqn|§} that to heat on average each 
element of disc gas up to a temperature of around T « 2000 
K once every viscous timescale requires about / ^ 0.08 of 
the available accretion energy. Thus we require that 



{q,) ^ f xq+, 
and hence that 

h/V, 



p^f 



1/2 



tth 



3.3 X 10" 



(34) 



(35) 



We have argued that each chondrule heating event 
should take place on a scale of ~ /i « 4 x 10^^ cm. If we 
divide the disc locally (say 2-4 AU) into A^box boxes, each 
of size ^ h? then 



Abox ~ %R^.2H/h^ ~ 1.4 X 10^ 



(36) 



Therefore the number of these boxes which is 'hot' at 
any one time is 



Ahot ~ P X iVbox ~ 4.6 X 10" 



(37) 



Since the orbital period at this radius is P = 1.6 x 10* s, each 
box is therefore heated for a fraction /hot ~ (h/Vs)/P ~ 
3.6 X 10"'^ of an orbit. Similarly the average time between 
these large flare events is given by fflare ~ -P/Ahot ~ 100 yr. 

Assuming that the reconnection rate for loops is roughly 
the same as the rate at which they are stretched by the disc 
shear flow, so that on average each loop reconnects once per 
orbit, then the number of reconnections per orbit which give 
rise to major flares, and so chondrule heating, is 



Arcqu ~ Ahot//hot ~ 13. 



(38) 



7.2.3 Loop strength distribution and relation to the 
dynamo 

We have argued that only about / ~ 0.08 (eqn[9]) of the en- 
ergy is released in heating events which are powerful enough 
to cause chondrule formation. If we make the simplifying 
assumption that magnetic loops are all approximately the 
same physical size (~ /i) but have varying magnetic field 
strength, then this implies that most loops have only weak 
field strengths. This ties in with the argument (eQn l20p that 
to obtain a viscosity parameter of a r; 0.01 we only require 
an average disc field of order {B) ~ 0.5 G. 

If all the loops are the same size, then the total number 
of loops would be 



Aioop, ~ Abox ~ 1.4 X 10^ 



(39) 



If, in a steady state, each loop reconnects once per dynamical 
time (~ orbital period) then over this time there are a total 



of A, 



loops 



1.4 X 10 reconnections, of which Ars 



13 



give enhanced heating, and so chondrule formation. These 
strong fiares provide a fraction / ~ 0.08 of the total power of 
the disc, and so the energy released in a weak reconnection 



event is 



A, 



rcqu oops 



1.3 x 10 times that released 



in a strong one. The frequency distribution f{E) for events 
which release an energy E (cf. Charbonneau et al., 2001) is 
usually assumed to take a power law form 



f{E) oc E-", s > 0. 



(40) 



Then our reasoning above implies that s « 2.3. This is simi- 
lar to what is found for the energy release in solar flares, for 
which s is generally in the range 1.4 ^ s < 2.6 (Charbonneau 
et al., 2001). Note that s — 2 implies equal energy release at 
all scales, so that s = 2.3 implies that energy release occurs 
mainly (but marginally so) at small scales. For their incom- 
pressible loop dynamo model, Baggaley et al. (2009b) find 
s « 3, so that small scales strongly dominate. 



8 CONCLUSIONS 

We have argued that the need for strong intermittent heat- 
ing during the formation of chondrules in the protosolar neb- 
ula places important constraints on the way gravitational 
energy is released in accretion discs. Following Desch & Con- 
noly (2002) we suggest that this requires shock heating, and 
propose that the shocks occur as magnetic loops in the disc 
dynamo are released by reconnection events. In contrast, 
current shearing-box simulations of MRI-driven accretion 
disc dynamos predict largely subsonic velocity fields, and so 
leave the formation of chondrules unexplained. 

We have proposed that accretion disc dynamos may be 
much more inhomogeneous, in space and in time, than cur- 
rent simulations indicate. It would not be surprising if the 
input physics to current numerical simulation was incom- 
plete, and we have briefly discussed reasons why this might 
be so. We speculate that a better understanding of accretion 
disc energy release may require consideration of more com- 
plicated physical processes than are currently implemented 
in numerical codes. 
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